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Abstract 

We study the classical and quantum transport processes on some finite net- 
works and model them by continuous-time random walks (CTRW) and 
continuous-time quantum walks (CTQW), respectively. We calculate the 
classical and quantum transition probabilities between two nodes of the 
network. We numerically show that there is a high probability to find the 
walker at the initial node for CTQWs on the underlying networks due to 
the interference phenomenon, even for long times. To get global informa- 
tion (independent of the starting node) about the transport efficiency, we 
average the return probability over all nodes of the network. We apply the 
decay rate and the asymptotic value of the average of the return probabil- 
ity to evaluate the transport efficiency. Our numerical results prove that 
the existence of the symmetry in the underlying networks makes quantum 
transport be less efficient than the classical one. In addition, we find that 
the increasing of the symmetry of these networks decreases the efficiency of 
quantum transport on them. 
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1 Introduction 



Quantum walk(QW) as a generalization of random walk(RW) is obtained 
by endowing the walker with quantum properties [HIS]- The QW has been 
largely based on two standard variants, the discrete time QW(DTQW) [3] 
and the continuous-time QW(CTQW) [1]. In recent years, the DTQWs have 
been investigated on trees [4J , random environment [5] , single and entangled 
particles El E]- The CTQWs have been studied on line [TOl [III IE] , 
cycle [131 E] > one-dimension regular network [15] , regular graphs [16] , odd 
graphs [T7], cayley tree [THl [T21 [2D] , hypercube [21], small- world network 
star graphs [23^ I24j. dendrimer [25) . restricted geometries [26j, Apollonian 
network f27] and one-dimensional and two-dimensional networks with pe- 
riodic boundary conditions [28^ I29j . Moreover, the decoherent QWs have 
been considered on hypercube [3^, cycle [311 [32], long-range interacting cy- 
cle \33\ [M] and one-dimension regular network [35] . 

Since the QWs may exploit quantum mechanical effects such as entangle- 
ment and interference, the only special experimental techniques can be can- 
didate to implement them. Recently, some experimental implementations 
of both QW variants have been reported e.g. on microwave cavities [36] . 
ground state atoms [37|, the orbital angular momentum of photons [38] . 
waveguide arrays [39] or Rydberg atoms [301 1^ . 

Since RWs generate the classical algorithms in computer science and model 
the diffusion phenomena and non-deterministic motion on the complex net- 
works [32], it can be expected that their quantum extensions(QWs) pro- 
vide tool to implement quantum computing and to model quantum pro- 
cesses [12] . In quantum computing, it is a primary goal to determine 
when quantum computers can solve problems faster than classical com- 
puters [j^ 145] . Grover in 1996 proved that a classical computer requires 
0{N) steps to finding marked item among N items, whereas quantum 
computer can solve it using only 0{^/N) steps [46j. Several implemen- 
tations to Grover's algorithm have been reported by CTQW [H Wl\ [38] 
and DTQW \4:9\ \50\ 151] . To implement Grover's algorithm, the all items 
must be accessible by local moves. For instance, if items are located on the 
one-dimensional line, traveling from one end of the line to the other end 
requires moves, thus the classical or quantum algorithms can not find 
a marked item in less time than 0{N) [47J. It is interesting to consider a 
one-dimensional line and then gradually to change its structure so that, in 
every step, items on it to be accessible by the lesser moves, as shown in Figs. 
l(a-e). 
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Figure 1: Figs. l(a)-(e) show some networks with ten nodes and nine edges. 

On the other hand, in recent years the complex networks have been ap- 
phed to model very diverse systems in nature such as lattices in solid state 
physics and condensed matter [521 [531 [Ml ES] . this modeling process, the 
components of system and interactions among them are considered as nodes 
and bonds of network, respectively. It is difficult to study the complex net- 
works having many nodes or bonds or both. But the most complex networks 
are constructed from simpler networks which the number of nodes and the 
number of bonds are smaller. Since these simpler networks can successfully 
explain some features of whole network, there is an upsurge of interest in 
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studying of their various aspects. Hence, by studying of RW and QW on the 
above mentioned simple networks, we also highlight the interplay between 
the diffusion process and the network symmetry. 

In this paper, we first calculate transition probabilities between two nodes 
of network for the CTRWs and CTQWs, then we numerically show that 
for CTQWs on the mentioned networks due to interference phenomenon, 
even for long times, there is high probability to find the walker at the initial 
node. Thus, the quantum transition probability distribution on underlying 
networks depends significantly on the initial node. To get a global informa- 
tion about efficiency of walk, we average the return probability over all nodes 
of network. Then, we evaluate the transport efficiency on networks(a-e) by 
the rate decay and by the asymptotic value of average return probability. 
The numerical results show that the existence of symmetry in the mentioned 
networks causes the quantum transport to be less efficient than the classical 
counterpart. In other words, the results denote that the increasing of the 
symmetry of these networks decreases the efficiency of quantum transport 
on them. 

Our paper is structured as follows: after a brief summary of the main con- 
cepts and of the formulae concerning CTQWs in Sec. 2, we study the quan- 
tum transition probabilities and their long time average on the mentioned 
networks, in Sec. 3. In Sec. 4, we especially focus on the average return 
probabilities and the efficiency of classical and quantum walks. Finally, in 
Sec. 5 the conclusions are presented. 



In general, every network can be characterized by a graph G that consists 
of a finite nonempty set V oi N vertices together with a prescribed set X 
of q unordered pairs of distinct vertices of V so that each pair x = {u, of 
vertices in X is a edge of G [56] . 

Algebraically, a graph can be described by its adjacency matrix A whose 
elements are A{i,j) = 1 for G X and A{i,j) = otherwise. The 

Laplacian operator is then defined as L = Z — A, where Z is a diagonal 
matrix and Zjj is the degree of vertex j. Classically, the evolution of CTRW 
is governed by the following master equation |57j : 



2 CTQWs on networks 



dt 



N 




1=1 
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where Pkj{t) is the conditional probabihty to find the walker at time t 
at node k when starting at node j. 

The matrix T is the transfer matrix of the walk, T = (Tkj), and relates to 
the Laplace matrix through T = — jL, where we assume an unbiased CTRW 
so that the transmission rates of all bonds are equal and set 7 = 1- 
The CTQW as the quantum mechanical extension of CTRW is obtained by 
replacing the Hamiltonian of the system by the classical transfer matrix, 

F = -T yy. 

We denote the state associated with node j of network as \j) and take the 
set {\j)} to be orthonormal. Thus, the solution of Eq. (1) is 

P,,{t) = {k\e'^\j). (2) 

Quantum mechanically, the states \j) span the whole accessible Hilbert 
space. The time evolution of state \j) starting at time to is given by 
\j,t) = U{t,tQ)\j), where to) = exp[—iH{t — to)] is the quantum me- 
chanical evolution operator. Hence, transition amplitude akj{t) from state 
|j) at the time to state \k) at time t is 

a,j{t) = {k\e-'^'\j). (3) 

From the Schrodinger Equation(SE), we obtain 

d ^ 
ij^^k,j{t) = J2HMaij{t), (4) 

1=1 

where h = 1. We assume that £"„ and \qn) denote the nth eigenvalue and 
eigenvector of H, respectively. The classical and quantum transition prob- 
abilities between two nodes can be written as |58j . 

N 

^fc,iW = ^e-*^"(fc|g„)(<?„|j) (5) 

n=l 



n=l 

Note that the normalization condition for Pkj{t) is 'Yl!k=i-^k,j{i) = 1 and 
for ak,j{t) is Y.k=i = 1- 

Classically, the transition probabilities converge to the equip-partitioned 
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probability jj^, whereas the quantum probabihties do not reach any finite 
value but after some time fluctuate about a constant value. This value is 
determined by the long time average(LTA) which is defined by using Eq. 

(6) m-- 



Xk,j 



1 f'^ 

lim - TTk,j{t)dt = y2^Er,,E,„{k\qn){qn\j){j\qm){qm\k), (7) 

>00 1 In 



where SE„,Em = 1 ^^v En = Em and 6E„,Em = otherwise. Since some 
eigenvalues of H may be degenerate, the sum in the Eq. (7) contains terms 
belonging to different eigenvectors. 



3 Quantum transition probabilities 

In this section, we illustrate the importance of choosing the initial node in 
CTQWs process on the networks. For this aim, we first consider the quan- 
tum transition probabilities in intermediate times when the walker started 
from a special node. Then, we generalize our result for long times and the 
whole nodes by studying the LTA probabilities. 

To obtain the quantum transition probabilities (see Eq. (6)), we need all the 
eigenvalues and eigenvectors of H, thus we make use of standard software 
MATLAB. 

In Fig. 2(a-e) we present the quantum transition probabilities TTk,j{i) for 
networks(a-e) where j = 5,4, 1, respectively. Figs. 2(a-e) show the high 
values for ir^^^, ir^^i, 713^3, 712,2 and vri^i. In other words, these figures prove 
that at short times (from t = Is to t = 20s), the quantum return probability 
is large and thus the CTQWs hold such dependence significantly on given 
starting node. 

Now it is natural to pose the following questions: 

Whether the such behavior continues at long times? 

Do the CTQWs hold such dependence on the other starting nodes? 

To address these questions, in Figs. 3(a-e) we plotted the LTA probabilities 

for all nodes k and j pertaining to networks(a-e), respectively. 

The axes x and y show nodes k and j of network respectively, and Xk,j is 

presented on the axis z. Figs. 3(a-e) show the symmetry characterizing the 

quantum transition probability, namely Xk,j = Xj,k- This can be derived 

directly from Eq. (7), recalling that H itself is symmetric and real. 
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In Figs. 3(a-e), X5,5, X4,4, X3,3, X2,2 and xi,i have the high peaks, respec- 
tively. These peaks are a consequence of the constructive interference due to 
reflections at peripheral sites and boundaries. They prove that, even at long 
times, the quantum probabilities depend on significantly to the starting node 
which is consistent with findings reported in [2 m i29 t [^ for Cayley tree(CT) 
and square torus(ST). Hence, to get global information(independent of the 
initial node) about CTQWs on networks(a-e), we must average over the 
whole network nodes. 



4 Efficiency of classical and quantum walks 

The transport processes implemented by CTRWs and CTQWs will be more 
efficient if the walker rapidly spreads over the network i.e. there exists a 
fast delocalization [26]. On the other hand, a fast delocalization implies 
a small probability to return(or stay) at the initial node, thus the return 
probability can be good candidate to evaluation the transport efficiency. 
But according the result of Sec. 3, to get a global information about the 
transport efficiency, we need a global quantity which is independent of the 
initial node. For this aim, we average the return probability over all nodes 
of network and obtain the average return probability. 

Hence, the classical and quantum average return probabilities are calculated 
by the following simple expressions, respectively [58] 

TV N 

P(') = NEP^At) = j^E'-''-' (8) 

j=l n=l 
N N 
j=l m,n,j=l 

It is evident that P{t) depends only on the eigenvalues and not to the 
eigenvectors of T while tt (t) depends on both the eigenvalues and the eigen- 
vectors of H. In the quantum case, we can define a lower bound for 7f(i) 
which depends only on the En, i.e,: 

1 ^ 

vf(t)>|aWp=|-^e-^"*P (10) 

n=l 

where a{t) = ^ J2f=i oijj{t). 

For the CTQW on a simple network with periodic boundary conditions, the 
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quantum average return probability equals to its lower bound i.e. 7r(t) = 
|Q;(t)p. The problem of the CTQW on hypercube lattices in higher d- 
dimensional spaces separates in every direction, thus we have a^'^j = [ajj{t)]'^ 
which results in Tr{t) = \a{t)\'^, too [59 j . 

We denote the degeneracy of eigenvalue En by Dn and rewrite Eqs. (8), 
(10) as 

En 

m > Ht)f = ^ E ^nI),ne-'(^"-^-)*. (12) 

According to Sec. 2, finally Tt{t) fluctuates about a stationary(asymptotic) 
value given by x- 

X= lim ^ r Tr{t)dt = ^ ^En,Ej{j\qn)Wj\qm)?. (13) 

n,m,j 

By Cauchy-Schwarz inequality or by taking the LTA of |Q;(t)p, we can obtain 
a lower bound of x which does not depend on the eigenvectors |60j : 

Xlb= ^YSEn,E„,- (14) 

n,m 

In fact this equation provides the exact asymptotic value of |a(t)P and the 
lower bound of asymptotic value of 7r{t). Since some eigenvalues of H might 
be degenerate, the above sum is equal to the number of non-degenerate 
eigenvalues plus the number squared of degenerate eigenvalues. 
The relation among the average return probability and the efficiency of 
transport can be considered from two different points of view. From the 
first point of view, the decay rate of the average return probability is pro- 
portional to the transport efficiency. The reason being that a quick decrease 
of the average return probability results -on average- in a quick increase 
of the probability for the walker to be at any other but the initial node, 
which provides a more efficient transport. From the second point of view, 
the asymptotic value of the average return probability has a inverse relation 
with the transport efficiency, because a large asymptotic value of the average 
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return probability implies a large probability to return at the initial node, 
which means inefficiency. 

In [60], authors studied the decay rate of P{t) and |Q;(t)|^ on the networks 
with two distinct eigenvalue spectrums: uniform degeneracy and one highly 
degeneracy. 

They numerically show that for the large networks (A^ S> 1) whose eigenval- 
ues have uniform degeneracy, the quantum walk can be more efficient than 
the classical one. Also, for the networks whose only eigenvalue Ei has high 
degeneracy Di, they used the following approximate equation 

^]^[A' + 2 Yl DnDiCOs{{En-Ei)t)]. (15) 

By this equation, they found that does not show a decay to values 

fluctuating about 1/N but rather to values fluctuating about 1 — thus 
the quantum transport is less efficient. 

In the following, we study the classical and quantum efficiency on the 
small networks(with few nodes) mentioned in Figs, l(a-e). Note that all 
networks(a-e) consist of 10 nodes and 9 links. The symmetry of network 
is given by whose nodes having the similar situations. For example, nodes 
8,9,10 in network(b) and nodes 5,6,7,8,9,10 in network(d) have the similar 
situations, thus network(d) is symmetrically higher than network(b). On the 
other hand, the numerical determination of the eigenvalues of networks(a-e) 
indicates that network(a) has no degenerate eigenvalue while networks(b-e) 
have one degenerate eigenvalue 1 with the order of degeneracy 2, 4, 6 and 
8, respectively. Hence, in networks(a-e) the increasing of network symme- 
try results in the increasing of degeneracy of eigenvalue 1. We denote the 
degeneracy of eigenvalue 1 by Di and represent it as the degree of network 
symmetry. We divide the problem into two separate cases as follows: 

4.1 non-degenerate eigenvalues 

Here, we consider network(a) with non-degenerate eigenvalue spectrums. 
For this network, Eq. (12) can be written as 



Prom the above equation, one can infer that after long time the only 
terms with E^ = E^ contribute to the sum and therefore, the equation 
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becomes of order 0{j^). 

Fig. 4(a) shows the temporal behavior of P{t), Tt{t) and for net- 

work(a). We can see that the classical curve(bluc line) does not show the 
constant value at the intermediate times whereas after t ^ N = 10s, not only 
|a(t)p(soUd red line) but also 7f(f)(green line) fluctuate about the equiparti- 
tion value 1/10. As mentioned above, the exact asymptotic value of |a(t)P 
and the asymptotic value of 7f(t) can be reproduced by Eq. (14). Since the 
eigenvalue spectrum of network(a) has no degenerate eigenvalue, the only 
the number of non-degenerate eigenvalues contributes in Eq. (14), resulting 
in the value 1/10. To determine the scaling behavior of P{t) and |Q(t)p, 
we use the dashed block(t~"'^/^) and red{t^^) lines, respectively. These lines 
show that the quantum return probability decreases faster than classical one, 
thus the quantum walk is more efficient than the classical random walk. 

4.2 degenerate eigenvalues 

In the following, we consider the transport efficiency on networks(b-e) with 
degenerate eigenvalue spectrums. Figs. 4(b-c) represent the temporal be- 
havior of P{t), |a(i)P and 7t{t) for networks (b-e), respectively. 

In all the figures, for all times, one can see that the quantum average 
return probability (green curve) is higher than the corresponding classical 
probability (blue cTirve), and thus the quantum transport is less efficient 
than the classical one. 

On the other hand, from Figs. 4(b-e), we find that after some time the strong 
maxima of Tt(t) on highly symmetry networks(d,e) can well be reproduced 
by the lower bound |a(t)p. Also, the increasing of symmetry causes the 
dashed block curve of Eq. (15), which is plotted for Ei = 1, to become 
very close to the full expression for |a(t)p. For instance, in Fig. 4(e) the 
positions of the exact curve of |a(t)P and its approximate equation almost 
coincide. It means that the only highly degenerate eigenvalues contribute 
to |a(i)p, and there arc only slight deviations due to the remaining ones. 
Therefore, for highly symmetrical networks(d,e), we can well apply Eq. (15) 
which proves that, always, |a(t)p fiuctuate about 1 — 1/A^ which is larger 
than the classical equipartition 1/A^. Thus, the classical transport is more 
efficient than the quantum one which is in agreement with the above result. 
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Now, to study the details of the effect of symmetry in the efficiency of 
classical and quantum transport, we study the asymptotic values of P{t) 
and 

Network(b) have eigenvalue 1 with degeneracy 2-fold, (i.e. Di = 2). In Fig. 
4(b), after t ~ 30, we observe that P{t) reaches equipartition 0.1, while 
|d!(t)p fluctuates about the asymptotic value 0.12. Since the eigenvalue 
spectrum of network(b) includes eigenvalue 1 with degeneracy 2 and eight 
non-degenerate eigenvalues, Eq. (14) gives the asymptotic value of |a(i)P 
as 0.12 which confirms the numerical result. 

In network(c), the degeneracy of eigenvalue 1 is 4-fold and the other eigen- 
values are non-degenerate. In Fig. 4(c), we see that at about t ~ 12s 
the classical curve reaches the equipartition probability 0.1 while the lower 
bound |a;(t)P fluctuates about the saturation value 0.22, as can be derived 
from Eq. (14). 

Network(d) have eigenvalue 1 with degeneracy 6 and other eigenvalues non- 
degenerate. In Fig. 4(d), one can see that the classical curve reaches 0.1 at 
t ~ 20s, whereas the lower bound |a(t)P fluctuates about the asymptotic 
value 0.4(as can be inferred of Eq. (14)). 

Network(e) can well be represented a graph star with 10 nodes whose eigen- 
values have three discrete values: eigenvalue 1 with degeneracy 8, non- 
degenerate eigenvalues and 10 [58] • Fig. 4(e) shows that after t ~ 8s, 
classical probabilities reach 0.1, while not only |a(t)P but also |vf(t)p fluc- 
tuate about a constant value 0.66(as can be obtained from Eq. (14)). 
We can conclude that, on symmetrical small networks such as those an- 
alyzed here, CTRWs can spread faster than their quantum counterparts. 
Moreover, the classical curves(blue lines) are flattened out and P{t) tends 
towards a limiting value(-^) such that with increasing the network symmetry 
this asymptotic domain is obtained more quickly, which implies a more effi- 
cient classical transport, except for networks(c),(d), i.e. however network(d) 
is more symmetrical than network(c)(D; of network(d)>D; of network(c)), 
the limiting value for it occurs more late. 

But with increasing the degree of network symmetry(Z)i), the quantum av- 
erage return probabilities fluctuate about a larger saturation value, which 
implies a less efficient quantum transport, and there is not any exception 
among networks (a-e). 

The reason of these different behaviors can infer from Eqs. (8,9). Based 
on Eq. (8), in the classical transport the eigenvalues of the transfer ma- 
trix itself dominates the average return probability and the degeneracy of 
eigenvalues directly dose not play role in increase or decreasing P{t), while 
quantum mechanically (Eq. (9)), the degeneracies of eigenvalues are governing 
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■7f(t). Therefore, as the above numerical results showed, while the efficiency 
of classical transport on the mentioned networks has not any exact relation 
with symmetry, increasing the degree of symmetry decreases the efficiency 
of quantum transport on them. 

5 Conclusions 

In summery, we studied CTRWs and CTQWs on some networks with few 
nodes. We calculated the classical and quantum transition probabilities on 
the networks and by numerical analysis found that there is high probabil- 
ity to find the walker at the initial node for the CTQW on the mentioned 
networks due to interference phenomenon, even at long times. Thus, to get 
information about the transport efficiency, we averaged the quantity of the 
return probability over the all nodes of network. Then, we studied the effi- 
ciency of transport on the mentioned networks by studying the decay rate 
and the asymptotic value of the average return probability. The numerical 
results showed that the existence of symmetry in the mentioned networks 
causes the quantum walks to be less efficient than their classical counter- 
parts. Moreover, we found that the increasing of the symmetry of these 
networks decreases the efficiency of quantum transport on them. 
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